Coherent Disintegration and Stability of Vortices in Trapped Bose Condensates 
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We consider the intrinsic stability of the vortex states of a pure Bose- Einstein condensate confined 
in a harmonic potential under the effects of coherent atom-atom interaction. We find that stable 
■ vortices can be supported, and that vortex stability can be controlled by changing the inter-particle 

interaction strength. At unstable regimes, a vortex will spontaneously disintegrate into states with 
different angular momenta even without external perturbations, with the lifetime determined by its 
imaginary excitation frequencies. 
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I. INTRODUCTION 

(— > 

a 

^ ■ Vortices and their motions have long been an important branch of fluid mechanics. With the discovery of superfiuid 
helium II, a new idea developed: that the circulation in a superfiuid vortex must be quantized jjj. The consequences 
of these quantized vortices is profound and the understanding of vortex dynamics plays a key role in the current 
understanding of superfluidity. Moreover, the detection of individual singly quantized vortices has vividly established 
the true macroscopic quantum nature of these remarkable degenerate fluids. Intimately related to the observation 
of the superfiuid state of 4 He is the evidence for the concurrent formation of a Bose-Einstein condensate (BEC) p). 
An important challenge is to clarify the link between these fundamental and important phenomena — Bose-Einstein 
£j ' condensation, superfluidity and the formation of a macroscopic quantum state. 

The recent observation of Bose-Einstein condensation in dilute alkali vapors [^|-^| has presented a striking new 
system for investigation; that of the dilute degenerate Bose gas. The alkali BECs differ fundamentally from the 
helium BEC in several crucial ways. BEC in both bulk liquid helium and the dilute helium "gas" are free systems(the 
^vq ■ "gas" BEC is created by introducing helium into a porous glass known as Vycor Q). By contrast, the alkali vapor 
BECs, although free of container walls (and/or the Vycor host), are created within the confines of a trapping potential. 
There is another major difference: in the trapped alkali condensates, samples can be prepared in which essentially all 
of the atoms are Bose condensed. By contrast, in bulk superfiuid 4 He, although the superfiuid fraction can be near 
unity, momentum distribution measurements have shown that the bulk condensate fraction is closer to 0.1 with the 
remainder of the particles in finite momentum states. As researchers improve their ability to create and manipulate 
these new trapped gaseous condensates, a series of important questions naturally arise: Does the gaseous BEC support 
superflow? Is it indeed a superfiuid? Are there stable vortices? This last question is the subject of this paper. 

The problem of vortex state excitations has been recently treated by others. Sinha @ investigated the low-lying 
Y , modes under Thomas-Fermi limit and Dodd et al. Q obtained the normal mode spectrum of a single quantized vortex 
state as a function of the number of condensate atoms for a BEC confined in a TOP trap ||. However, the important 
question of vortex stability was not addressed by these authors. More recently, Rokhsar studied the stability 
properties of the trapped vortices and argued that vortices are unstable due to the existence of a bound state inside 
the vortex core. However, throughout his analysis, the transition from a vortex state to a core state requires the 
presence of thermal atoms which serve as a reservoir to conserve the energy and angular momentum in the process. 
Hence, the more fundamental question concerning the intrinsic stability of an isolated vortex (i.e., without the external 
influence of thermal atoms) remains unanswered. 

In this paper, we approach the problem by assuming that all atoms are in the condensate such that scattering with 
thermal background atoms can be neglected. This allows us to focus on the intrinsic coupling within and between 
different vortex states and on the effect of this coupling on vortex stability. (Here, we use the word "intrinsic" 
to emphasize coherent coupling between the condensate atoms.) We find that stable vortex states can in fact be 
supported and show that whether a vortex state is stable or not is determined by its angular momentum and the 
nonlinear inter-particle interaction strength. Furthermore, we point out that the lifetime of an unstable vortex can 
be directly determined from the frequencies of the collective excitations. 

The paper is organized as follows. In Section II, we discribe our physical model and define the stability criterion. 
Our main results are presented in Section III, where the stable and unstable regions of trapped vortices are identified. 
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We also present a physical interpretation of the meaning of the instability. Finally, we give a summary and compare 
our work with others in Section IV. 



II. PHYSICAL MODEL 



To simplify our calculations, we consider a condensate confined in a 2d isotropic harmonic potential with trap 
frequency coq at zero temperature. In current experiments, condensates are achieved in 3d traps with cylindrical 
symmetry. A quasi-2d situation can be realized when u>± <C lo z , where u± and w z are transverse and longitudinal trap 
frequencies respectively fic| . In this limit, one can produce a "pancake" -shaped condensate with all the atoms lying 
in the lowest harmonic oscillator state in z-direction and hence, the degree of freedom in z-coordinate is frozen. Our 
treatment lies within the Hartree-Fock-Bogoliubov approximation. First, we calculate the macroscopic wavefunctions 
of the condensate in a vortex state. Next, we find the collective excitation frequency u) of the state. The stability 
criterion is defined as lm(cj)=0 IP]] , selected because if Im(u>) 7^ 0, then fluctuations around the condensate can grow 
exponentially in time and hence induce instability. 

In the Gross-Pitaevskii treatment [^2) , the energy for N condensed bosons of mass m is given by the functional: 

ffl^ = J dv (**f * K + V |* K | 2 + \NU\^ K f) , (1) 

Here 

tf„(r) = t> K (r)e lfi6 , k = 0, ±1, ±2, ... (2) 

represents the wavefunction of the macroscopic vortex state with azimuthal angular momentum nTi. T = 
— h 2 \7 2 /2m, V = mu^r 2 /2 are the kinetic and potential energy operators, respectively, and the coupling constant 
U describes the interactions between condensate atoms. In the quasi-2d situation considered here, the coupling con- 
stant takes the form U — ^^/ttTiuj z ^ z a jl3), where £ z = y / h/2muj z is the harmonic oscillator length in z-dimension. In 
our analysis, the solution that minimizes Eq. (1) is found iteratively using a finite elements method (FEM) p^ , p^| . 
In our calculations, we normally used 20 elements, with 2 nodes and 3 degrees of freedom for each element. This 
numerical method is very efficient, and it typically took no more than a few minutes to find the wavefunction \1/ K in 
a Cray-YMP2E/232. 

With the solution of fy K at our disposal, we can now calculate collective excitation frequencies by solving Bogoliubov 
equations B: 

(£ - hio x ~ fx K ) u x (r) + NU[* K (r)} 2 v* x (t) = 0, (3a) 
NU[^l(r)] 2 u A (r) + (£ + huj x - (m k ) v* x (r) = 0, (3b) 

where (j. K is the chemical potential for state ^ K (r), u x (r),v x (r) are normal mode functions with mode frequency lu x , and 
C = T + V + 2NU\^ K (r)\ 2 . It is straightforward to show that if \E , K (r) is given by Eq. (2), then u x (r),v x (r) must have 
definite angular momentum compositions k u % and n v Ti respectively such that u x (r) = u x {r)e lli ^ 9 , v x (r) = v x (r)e lKv6 , 
and k u + k v = 2k S. 



III. RESULTS AND INTERPRETATION 



Eqs. (3) were transformed to an eigenvalue problem for a finite-sized matrix and solved using the FEM. Our goal 
here is to find mode frequencies with non-zero imaginary part, in order to determine the vortex stability. As in 
the case of ground state, the vortex stability properties for a condensate with repulsive inter-particle interaction are 
drastically different from that for a condensate with attractive interaction. We will discuss these two cases separately. 

Repulsive interaction, i.e., U > — When we calculate the collective excitation frequencies of a single quantized 
(re = 1) vortex state $1, we find that all the excitation frequencies are real, which means that ^1 is always stable. 
Next, we consider a double quantized (k = 2) vortex state ^>2- Here we find that complex frequencies only exist for 
k u = and k v = 4 (Without loss of generality, we assume k > 0, and k v > k > k u .). We find that for any other pairs 
of (k u , k v ), the excitation frequencies are all real. Furthermore, for values of NU for which the vortex is unstable, we 
find that there exists at most one complex frequency. Fig. 1(a) shows the imaginary part of the complex frequency 
as a function of interaction strength NU. As we can see, in this particular channel [i.e. choice of (k u , k v )], the 
parameter space of NU is divided into alternating stable and unstable regions. In the unstable regions, the inverse of 
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lm(u}) determines the lifetime of the unstable vortex. For the parameter range described in Fig. 1, the most unstable 
vortex state will decay after several periods of the harmonic trapping potential. We stress that the details of how 
the condensate will evolve under these instabilities is beyond the capability of the mean-field treatment and requires 
further investigation. 
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FIG. 1. (a) Imaginary part of the complex frequency of a double quantized vortex state ^2 as a function of interaction 
strength NU, for k u — and k v — 4. (b) Same as (a) for a triple quantized vortex state $3. Solid line: k u — and n v = 6; 
dashed line: k u — 1 and K v = 5. Frequency is in units of trap frequency loq, U in units of (Hluo^ 2 ), where £ = (h/2mujo) 1 '' 2 is 
the harmonic oscillator length. 



For a general state ^f K , our numerical calculations show that there are (k — 1) unstable channels that possess 
complex excitation frequencies; those with k u — 0, 1, ...,«; — 2 and k v — 2k — n u . Fig. 1(b) shows the imaginary part 
of the complex frequency for a triple quantized vortex state ^3. We can see a similar pattern as in Fig. 1(a), but 
here, there are two unstable channels. Each channel shows its own quasi-periodic behavior as a function of NU. The 
two channels have quite different "period" and characteristic width of unstable regions. At first look, this may appear 
rather unexpected. To interpret this behavior we will show that each unstable region in NU-sp&ce represents a decay 
channel in which two atoms from the given vortex state scatter into two new states, with angular momenta K u h and 
K v h, respectively, thus inducing instability in that initial vortex state 

First, let us define a boson field operator as: *f(r) = \/]V , I' K (r) + $(r), where the c-number ^ K (r) denotes the 
one-body wavefunction for the condensate and ?/>(r) is the field operator for the fluctuation part The second 

quantized Bogoliubov Hamiltonian reads: 

K B = J dv ^{v)[C- ^{v) + [^NU J drft(j.)ft(r)* K * K + h.c.}, 

where the c-number part independent of ip(r) has been neglected. We can further decompose ip(r) as ip(r) = 

Y] a n ,a^ii,Q(r), where a„ !( , is an annihilation operator associated with a single-particle state <j) n ,a- The set of states 
{4>n.a} is defined as the eigenvectors of C with eigenvalues e n ,a, he., C4>n,a = £n,a<fin,a, with subscripts (n, a) labeling 

the radial and angular quantum number respectively. Hamiltonian Kb may then be rewritten as: Kb = Hq + Hi, 
where 

H 

and 



= H A ( n «'««; n «'««) a l u ,«„ a L,K« + h -' 



n u ,K u n v .k v 
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A(n u ,K u ;n v ,K v ) = -NU J dr 0* u ^r)0*^(r)# K (r)*, t (r) . (4) 
In the interaction picture, a^ n a (t) — a\ L a e , ( e "' o_ ' 1 ' i ) | ] and the Hamiltonian is given by: 

E E A K, «„)e i(e "--"- +e "»-»- 2 ' 1 -M„,«„4„,K„+'i-c. (5) 

A (n u , K u ]n v , n v ) is non-zero only when k u + k„ = 2k, which is a direct consequence of the conservation of angular 
momentum. The interaction described by Hamiltonian (|^) is analogous to parametric processes in quantum optics 
where instability can occur under certain conditions. For example, one can build up large numbers of photon pairs 
(signal and idle) from the vacuum via parametric down conversion if the field frequencies satisfy a parametric resonance 
condition. In our case, the fluctuation in mode pair (0n u ,re„) 4>n v .n v ) grows exponentially when 

l e n„,«„ + - 2^ K | < A (n u , k u ; n v , k v ) , (6) 

and hence the vortex is unstable under such resonance condition. We emphasize that the instability implied in this 
picture is purely quantum mechanical. The atoms in the vortex can spontaneously disintegrate into 0n u . K „ and ni ,. Kt , 
states without the need of external (classical) perturbations, such as the interaction with the thermal background 
gases or perturbation of the trap. 

For a k = 1 vortex, our numerical calculations show that there exists no particle states that satisfy the resonance 
condition (||), in support of our prediction that a single quantized vortex state is always stable for U > 0. For a 
k = 2 vortex, we find a pair of particle states (0o,o 7 00,4) indeed satisfy inequality (^|). In the weak coupling limit, 
we can calculate the decay rate of the double quantized vortex state $2 using the Hamiltonian (||) by neglecting all 
the nonresonant terms (i.e., keeping only terms with n u = n v — 0, k u = 0, k v = 4 ). The results are shown in Fig. 2 
along with the imaginary part of the complex excitation frequencies of vortex state We can see a clear qualitative 
agreement between the two results. The agreement can be significantly improved if the contribution from states 
(</>i,0j 0i, 4) is also included in calculating the decay rate (see Fig. 2). We remark that although useful for interpreting 
our vortex stability predictions, the parametric resonance picture is valid only for the weak interaction regime. This 
is because a strong interaction can drastically change the frequencies of the oscillators and introduce mixing among 
different particle states. Further work would be necessary in order to understand all aspects of the complex structure 
shown in Fig. 1, particularly for large NU. 
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FIG. 2. Solid line: imaginary part of the complex frequency of a double quantized vortex state ^2 as a function of interaction 
strength NU for k u — and k v — 4; dashed line: decay rate of the vortex state ^2 as a function of NU calculated using 
Hamiltonian (4), only contributions from the resonant states (0o, 0,^0,4) are included; dot dashed line: same as dashed line, 
but include contributions from two more states (01,0,01,4). 

Attractive interaction, i.e., U < — A condensate with strong attractive inter-particle interaction is known to be 
subject to collapse. However, a metastable condensate with a small number of atoms can still exist 0,0- Fig. 3 shows 
the imaginary part of the complex excitation frequency for a single and a double quantized vortex state as functions of 
NU. Fig. 3(a) shows that ^1 is stable for sufficiently small attractive interaction, but unstable for larger interaction 
strength. For ^2, as we can see from Fig. 3(b), the channel (n u — 0, k v — 4) possesses complex frequency for all 
negative values of NU instead of showing a quasi-periodic pattern as in the case of repulsive interaction. Furthermore, 
we find that, similar to ^1, other channels which are stable for NU > become consistently unstable for sufficiently 
large \NU\ [We only show two such channels in Fig. 3(b).]. Our calculations show that, for U < 0, stable vortices only 
exist for a single quantized vortex state in the weak interaction regime [see Fig. 3(a)]; multiple quantized vortex state 
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(i.e., k > 1) is always unstable. It has been speculated that the existence of vortices may help stabilize a condensate 
with negative scattering length flgf] . However, as we show here, although such vortices may seem to be more stable 
against the collapse when compared to the ground state, they remain fundamentally unstable and small fluctuations 
will eventually destroy such vortices. 



s 



\ \ 

\ 1 \ 
\ \ \ 




(a) 


v» \ 




(Ku, Ku) 


- V \ 




(-2,4) _ 






(-1 , 3) 






- - (0, 2) 


\ 

\ 

\ 







-50 -40 -30 



10 




-100 -80 -60 -40 

NU 

FIG. 3. Imaginary part of the complex frequency of a single (a) and a double (b) quantized vortex state as a function of 
interaction strength NU, with NU < 0. 



IV. SUMMARY AND DISCUSSION 



In summary, we have calculated the collective excitation frequencies of a Bose-Einstein condensate in a vortex 
state and have established intrinsic stability regions for these vortices. We have shown that, even without any 
perturbation, an unstable vortex can still decay spontaneously. For repulsive inter-particle interaction, we found that 
single quantized vortices are always stable, while imaginary excitation modes divide the interaction energy axis (NU) 
of multiple quantized vortices (k > 1) into alternating stable and unstable regions. Hence, one can control the vortex 
stability by varying the value of interaction strength, which in turn can be achieved by changing the scattering length 
p9| , p0[ , particle number or trap frequency. This provides us with the possibility of studying condensate evolution 
under the effect of imaginary modes. 

For a condensate in vortex state, there may exist quasiparticle states with negative frequencies. One such negative 
frequency state was identified by Dodd et al. in Ref. || . The presence of negative frequencies implies that there exist 
states with lower energy. However, this does not necessarily mean that the condensate is unstable if no mechanism 
exists to drive the system to these lower energy states pl)| . In Ref. ||, Rokhsar considered the instability arising 
from the incoherent interactions between condensate and thermal atoms, which induce the transition to the negative 
frequency core state. In contrast, in the present paper, we study the intrinsic stability of vortices in a pure condensate, 
excluding such incoherent processes while focusing on the coherent interactions within the condensate. In our work, 
instability occurs as a coherent process such that an unstable vortex state will disintegrate into different angular 
momenta states. We found that stable vortices can be supported in harmonic traps as long as the temperature is 
low enough such that the effects of thermal atoms are insignificant. At temperatures when thermal atoms cannot be 
neglected, both coherent and incoherent processes will be present and each will have its effect on vortex stability. It 
remains to be seen which process will be dominant. Further investigations should also include the possible influence 
of trap anisotropy and the dynamics of the disintegration processes. 

Recently, vortex stability in 2d harmonic trap is studied by Caradoc-Davies et al. through a direct numerical 
simulation p2fl . In that study, a blue detuned laser beam is applied to perturb the condensate in a vortex state. They 
found that the single quantized vortex is indeed stable, while a double quantized vortex can disintegrate into unit 
vortices under external perturbation. These results are consistent with ours presented in this paper. 
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Finally, as an example, let us consider a 23 Na condensate (scattering length a w 3 nm) in a harmonic trap with 
w_l = 2?r x 10 Hz and w z = 2n x 200 Hz, in units of fiui^, U w 0.02. The plotted range of W from to 4000 in 
Fig. 1 corresponds to particle number ranging from to 2 x 10 5 , well within the capability of current experiments. 
Recently, several methods on how to generate vortex states in alkali atomic BECs have been proposed [j23| - With 
current technology and fast progress on this field, our study on vortex stability should be experimentally testable in 
the near future. 
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